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Abstract 

Using the QCD dipole picture of the hard BFKL pomeron, we derive the general 
expressions of the elastic and inelastic components of the proton diffractive structure 
functions. We obtain a good 7 parameter fit (including a secondary reggeon contribution) 
to data taken at HERA by the HI and ZEUS collaborations. The main characteristic 
features of the model in reproducing the data are discussed, namely the effective pomeron 
intercept, the scaling violations and the beta dependence. A difference obtained in the 
separate HI and ZEUS fits leads us to analyse directly the differences between both 
measurements. Predictions on R, the ratio of longitudinal to transverse photon cross 
sections are performed and lead to very large values at high beta and large virtuality Q 
which may lead to a discrimination between models. 



1 Introduction 



The HERA data for deep-inelastic scattering at high energy and high Q 2 contain a sizeable 
fraction of events in which one observes a large rapidity gap in the forward region J], |2j. 
These events result from a colour- singlet exchange between the diffractively dissociated virtual 
photon and the proton. Since this phenomenon is present even for high virtuality of the photon, 
it is called hard diffraction, at variance with soft diffractive phenomena. Several theoretical 
formulations have been proposed 0]. Among the most popular models, the one based on a 
pointlike structure of the pomeron [|J has been studied quantitatively using a non-perturbative 
input supplemented by a DGLAP evolution M. In this formulation, it is assumed that the 
exchanged object, the pomeron, is a colour-singlet quasi-particle whose structure is probed in 
deep-inelastic scattering. 

There exists a different approach |5|, || in which the cross sections are determined by the 
interaction between colour dipole states describing the photon and the proton. Indeed, 
it is well-known that the photon can be analyzed in terms of qq configurations || while it 
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has been shown || that the small-x structure function of the proton can be described by a 
collection of primordial dipoles with subsequent perturbative QCD evolution. More specifically 
H, the combination of the dipole description of perturbative QCD at high energy and the 
Good- Walker mechanism [T(| leads to a unified description of the proton total and diffractive 



structure functions |TTJ . However, a quantitative description of the data following || was still 
lacking. The aim of the present paper is to provide a simple and quantitative formulation fitting 
the experimental data. 

In the dipole approach, two components are shown to contribute to the diffractive structure 
function (see fig.l). First, an elastic component (fig.l-a) corresponds to the elastic interaction 
of two dipole configurations. It is expected to be dominant in the finite (3 region, i.e. for small 
relative masses of the diffractive system. Second, there is an inelastic component (fig. 1-b) where 
the initial photon dipole configuration is diffractively dissociated in multi-dipole states by the 
target. This process is expected to be important at small f3 (large masses). Note that both 
components are obtained in ref.|J through perturbative QCD resummation corresponding to 
the BFKL approach ||12|| . In these respects, it is based on a "hard pomeron" approach at 
variance with e.g. a non perturbative pomeron input evolved using DGLAP evolution |]]. 

We present here fits of the published diffractive data (T|, |2| based on a simple with 7 free 
parameters and analytic parametrization for the diffractive structure functions in the QCD 
dipole model. Both fits are successful showing that this approach is a good candidate for an 
understanding of hard diffraction. The main input of the calculation is a formulation of the 
elastic QCD amplitude T(r, p, b; xp) for the diffusion of a dipole of size r on a dipole of size 
p with rapidity gap log l/xp at fixed impact parameter b. This 6-dependent input allows one 
to factorize the different integrals involved in the perturbative calculation of the diffractive 
components. The dipole model for the proton || is used to provide predictions for the HERA 
data. We also include a phenomenological secondary Regge trajectory which is known to play 
a role in the limited domain of large mass and small rapidity gap ]T|| . 

The outline of the paper is as follows: in section 2, we introduce the dipole formalism and 
in particular the input elastic dipole-dipole amplitude T(r, p, b; xp) following the requirements 
of the BFKL dynamics. Using this amplitude, we derive the general expressions for the elastic 
(2.9, 2.10), and inelastic (2.15) diffractive components for both transverse and longitudinal 
polarizations of the photon. In section 3 we obtain an analytical form of the proton diffractive 
amplitudes used to fit the HERA data. In section 4, the resulting fit is displayed. The simple 
form of the amplitudes makes transparent the discussion of the model properties e.g. the relative 
contribution of the two components, the effective exponent for the rapidity gap dependence, 
the origin of the scaling violations, and the j3 dependence. Conclusions and outlook are given 
in the final section 5. 



2 The QCD dipole formalism 

As mentionned in the introduction, the main input of our approach to the structure functions 
is the interaction amplitude T(r, p,b; xp>) for a dipole of size p with a dipole of size r, at 
impact parameter b where logl/xp is the rapidity gap. The knowledge of the amplitude as a 
function of b is essential for the calculation of the two components of hard diffraction |J. It 
is thus compulsory to start with a correct BFKL amplitude in the whole impact parameter 
space which goes beyond the commonly used expression for the forward elastic amplitude 
rpiforward) = j d 2 b t (r, p, b; x P ) . For this sake, one uses the conformal invariance of the BFKL 
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kernel in transverse coordinates We shall propose such a 6-dependent amplitude by asking 
for the two following requirements: 

(i) The integral over impact parameter space boils down to the known exact expression for the 



forward amplitude T^ orward ^ pL 03 



(ii) It has the correct high impact parameter approximation Jf] including the modification by 



a scale factor which has been deduced from conformal invariance p~3 |. 

Using the inverse-Mellin transform which appears in the solution of the BFKL equations 
T"2"|j, the amplitude reads: 



T(r,p,b;x ¥ ) 



7 2^(1-27)^7)^— J % , (2.1, 



2 b A 

where the factor 16 is the scale factor determined by conformal invariance. 

a N 

A P ( 7 ) = — x( 7 ); x{l) = 2^(1) - V>(7) - ^(1 - 7) (2.2) 

7T 



is the 7-dependent QCD pomeron intercept |]l2"fl , and 

= i^t^f < 2 - 3) 

is the two-gluon exchange elementary dipole-dipole cross section. As already mentioned, the 
amplitude meets the abovementionned requirements, since: 

(i) the integration over b leads to the known expression for the forward amplitude (for r > p): 

,+oc ,+oo^_ k ^ T ( fe ) = g 2 f ^Lrf( 7 ) ( I^j 7 X " Ap W = T (forward) _ 
Jr/4 A/4 7 2?7T 

Note that the lower integration bound r/4 on b is required to recover the correct expression 
for the forward amplitude. In the absence of an exact conformal invariant solution for all 



impact parameter [14), this cutoff will play an important role in the determination of the hard 
diffraction components. 

(ii) A saddle point approximation of the integral over 7 in formula (|2.1|) gives the approximate 
expressions in the large b approximation used in ref.|| 0], up to the scale factor 16. This factor 
is strictly speaking not present in the large b approximation, but appears when one takes into 



account the global conformal invariance of the BFKL kernel This scale factor leads to 

a much more central impact parameter distribution than the approximate expressions used in 
ref.H 0, and plays a major role in diffraction, as we shall see in the following. 

Using formula ( |2.1| ) as an input leads to an analytic formulation of the diffractive amplitudes. 
In this section, we shall first consider hard diffraction on a dipole of given size r. It will be 
shown in the next section how, starting from this case, it is possible to derive phenomenological 
expressions for hard diffraction on a proton. One has the following identities for the diffractive 
structure functions [BJ: 

f d (r> 2 M 2, _ Q A 1 da T , L _ Q 2 1 da TiL 

t< TjL {Q ,Xp,M ) - A7T2(Xem(3x¥ dM2 ~ A7l 2 aemX / df3 > W 

where dux (resp. dai) is the differential cross section for the scattering of a virtual photon 
with transverse (resp. longitudinal) polarization, on a dipole. 
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a) Elastic component 

Using the QCD dipole model for the elastic component (see fig. l.a), one writes 

r (eQ 



daZ£ N c a em e 2 ^ 2 l2l f 1 



dM 2 



2tt 



Q 2 r d 2 h I dzf T , L (z)\T TjL (b;x r )t 

Jr/A JO 



where 



and 



T TtL (b;x P )= / dp pT(r,p,b;x P )K lfl (Qp)J lfi (Mp) 
Jo 



z 2 a-z) 2 (z 2 



(2.6) 



(2.7) 



Mz) 

f L (z) = iz 3 (l 



(2- 



The notation " X " is defined by X = (z(l - z)f/ 2 X. 

Following ref.||, the formulae ( |2.6| -2.8) describe the interaction of a qq configuration of 
transverse size p of the virtual photon with a dipole of size r. The input elastic amplitude is 
taken from ( |2.1| ). The photon wave function in terms of qq configurations is projected on its 
transverse (T) and longitudinal (L) polarization states which differ only by the factors /t,l(^) 
and by the order of the Bessel functions (TTSJ involved: J\ are associated to the transverse 
component, and Kq, Jq to the longitudinal one. z (resp. (1 — z)) is the longitudinal momentum 
fraction carried by the antiquark (resp. quark). Note that interference effects are present |j| 
which explain the product of Bessel functions of different kinds (K and J). The cutoff in impact 
parameter space appears twice in formulae ( |2.6|) and in ( p.7|) , first by the upper bound 46 on p 
and second by a lower bound r/4 on the impact parameter b. In practice, we will release the 
first bound, which will greatly simplify the calculations and allow to find an analytical form for 
the diffractive amplitude. This approximation is expected to change only the normalisations. 

Inserting expression ( [2.1| ) for the dipole-dipole amplitude in eq. ( |2.6| , |2.7| ), one obtains the 
results for F® L ( |2.5 ) involving integrals over the impact parameter b, the momentum fraction 
z and the gg-pair size p. Interestingly enough, these integrations can be factorized and exactly 
computed using in particular the change of variable p —> p = (z(l — z)) x l 2 p. After a tedious but 



straightforward calculation using well-known identities [15], the successive integrations give 



£^-«(|)Vn 



Q 



<ta 

2ztt 1 V 2Q 0V ^ / 



2-u 



~A P ( 7j ) 



x [] {(l-2 7 ,)d(7 J )r(2-7 J )r(3-7 i ) 2 F 1 (2-7 J ,-l+ 7j ;2;l-/3)} 

J'=l,2 

7i+72 



X 



(7i+72-l)(3-2( 7l + 72 )) 



5(7i+72, 1/2) , (2.9) 



p D{el) 



16 3 N r e 2 



Xp 7T Z 



33 if 



n 

J'=l,2 



dji ( Q ) 
2i7r \2Q ^p) 



2u 



X, 



-A P (7j) 



x n {(l-27,M7i)r 2 (2- 7j ) 2^(2-7,, -i+ 7j ;i; i-/5)} 

3=1,2 

1 



X 



3-2( 7l + 72 ) 



S( 7l + 72 ,l/2) , (2.10) 
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where we have introduced the scale Qo = 2/r. Note that these expressions involve two coupled 
inverse- Mellin transforms in 71 and 72, each corresponding to one elastic dipole interaction, 
(see fig.l-a). 

b) Inelastic component 

As shown in fig. l.b, the inelastic component stems from the following process: the initial 
qq state of the virtual photon develops a set of colour dipoles through cascading @ and the 
diffractive component is due to the interactions of two of the produced dipoles with a target 
dipole, each of these being described by the amplitude (|2.1|) . The cross section can be expressed 
as HI: 



dp s c Q 2 Jr/A b 2 \ ' J 2in \ 2 

where 



(2.11) 



i6 



dg 2 T(r, g,b;x ¥ ) 



2\ / "'K 2"*- '!C1 "IT /o in\ 

e ) = /_ — Q . 2 ( 2 - 12 ) 



£ Arrajg 



and 



# M = YkH r(3-7)r 3 (2-7)r(2+7)r(i+ 7 ) 

T[1> 167 2 (l-7) 4 r(4-2 7 )r(2 + 2 7 ) 

tt ( \_ vin) r(3-7)r 3 (2-7)r(2+7)r(i+ 7 ) 

8 7 (l-7) 3 (l+7)(2-7) r(4-2 7 )r(2+2 7 ) 1 " J 

with 

K( 7 ) = 3 F 2 (l- 7 , 1-7; 1,3/2; 1). (2.14) 

The dipole-dipole interaction is represented by formula (2.12) which can be interpreted as 
the average squared transverse sizes of the dipoles seen at the impact parameter b. Indeed, the 
quantity T(r, g, b; xp)/47ia 2 g 2 is nothing else than the density of dipoles of size g emitted from 
a dipole of size r at impact parameter b |7|, [13[ . The inverse Mellin-transform and its integrand 



factor V(7) (|2.14|) comes from the calculation of the 1 — > 2 dipole branching at the virtual 
photon vertex ||. The different prefactors appearing in and in formula (2.13) come 
from the transverse and longitudinal wave-functions of the virtual photon, respectively. 

Inserting the expression (2.12) into eq. ( j2.11 ) and integrating over the impact parameter b, 
one gets: 

D(m) _ 2 iv 2 / n \ dlj l ~ 2l U h )x- Ar ^ ] } 



\f 2 [ 2m 2 — 27j 
dl ( Q 



3= 

27 



x 



^(^(i-^Sfl- (7 ) . (2.15) 



2m \4QoJ ' 2-7-71-72 



The formulae (2.9, 2.10, 2.15) condense in an analytical form the perturbative QCD predictions 
for the diffractive structure functions of a dipole of size 2/Qq. 
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3 Proton diffractive structure functions in the dipole model 



In the previous section, we have obtained the expressions for the diffractive structure functions, 
when the target is a single dipole of size r = 2/Qq. Our aim in this section is to use these 
formulae to obtain the diffractive structure functions of the proton. We have thus to implement 
a model for the proton which could describe it as a collection of primordial dipoles. 

a) Elastic component 

Applying the steepest-descent method, the saddle-point is determined by the usual BFKL-type 
integrand (between the square brackets in (2.9, 2.10)). It reads: 

7 sp =i(l-a(x P )log^=) . (3.16) 

where a(x) = ir/(7((3)a s N c logl/x) is the /cj_-diffusion coefficient coming from the BFKL 
evolution |12|, [L6|]. We have to use a special treatment for the pole at 71 + 72 = 1 in formula 
(2.9), since it contributes for the asymptotic saddle points at 7! = 72 = 1/2 when x P — > 0. All 
the other singularities are not relevant, since they are far from the saddle points. 

Noting that in the kinematical domain of interest, 1 — 27 sp = a(xp) \ogQ/2Q y / ~j3 <C 1, and 
working at lowest order, the results are: 



Ft 12 x w a W log — t= e 2C W/3 

7r 2g vP 



x(3(l-(3) 



(3.17) 



7T 



xp 2 



(3.18) 



with 



= l + A P ('i) = l + ^£41og2. (3.19) 



The pole at 71+72=1 in results in an extra logQ/2Q V7^ which may be large. 

The formulae (3.17, 3.18) give an explicit expression for the elastic components of diffraction 
in the QCD dipole model. The xp- dependent prefactor (xp) plays the role of the (hard) 

pomeron flux factor already discussed in references || |TIJ. Revealing new features of the model, 
the present explicit calculation shows that the relevant scale factor responsible for the scaling 
violations is 2Q ^/J3 where Q is the scale of the proton-dipole amplitude. The /^-dependence 
is found to be non polynomial. 

b) Inelastic component 

Using again the steepest-descent method for evaluating the three integrals in 7, 71 and 72 of 
formula (2.15), we find saddle-points at 

ls P = ^(1 - a(P) log — ) , 7 J>P = ^(1 - a(xp)) . (3.20) 
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Note that there is no coupling of the integrals at the saddle points; however, one should take 
into account the prefactor zeros at jj = 1/2, see eq. (2.15). All in all, the result reads: 

,,,„■„, [2 fl\ N c e 2 a 5 s _2op+i 3/ x Q -^log 2 ^ i/m/P-Ap 



V 7T V2/ 7T 4 Q (3 

Comparing our expression ( |3.21| ) with the original result (formula (11) in the second paper of 
ref.0), we find identical results but for the scale of the transverse momentum in the BFKL 
evolution: 4Qo instead of Qq. This change comes from the scale factor determined by conformal 
invariance (see 2.1) and was not considered in the previous approximation II. This modified 
scale factor has an important phenomeno logical effect. 

In the QCD dipole model, there is a tight connection between the total and diffractive 



structure functions. In refs.P, 11 1, the QCD dipole model has been successfully applied to the 



proton total structure function in the small-x kinematical domain. A model for the primordial 
dipole configurations in the proton has been introduced through a distribution with average 
size and the non perturbative features of the proton target manifest themselves only through 
this scale and the global normalization constant related to the density of primordial dipoles (for 
a more complete discussion, see ref.||). We thus use this representation of the nucleon target 
in formulae (3.17, 3.18, 3.21 ) to derive the model for the proton diffractive structure functions. 
In this context, the scale parameter is reinterpreted as Qq = 2/ <r> and we introduce the 
arbitrary normalizations N? . , N^ l \ for F®l an d F^^ multiplying formulae (3.17, 

3.18, |3.21| ) respectively. These unknown normalizations reflect the non-perturbative primordial 
dipole configurations contributing to the various proton structure functions. 

With this model, we are ready to write a full parametrization adequate for the description 
of the data. The free parameters of the dipole model are o?p, which is related to the fixed 
coupling constant a s in the BFKL scheme at leading order (see eq. ( |3.19| )), Q , corresponding 
to a non-perturbative scale for the proton, and the three normalizations N!f , N^ l \ N^ n K As 
is now well-known, a secondary trajectory based on reggeon exchange is added in order to take 
into account the large-mass and small rapidity gap domain. Reggeon exchange can here be 
simply parametrized in the following way: 

F 2 D{R) (x P ,P,Q 2 ) = f(x F )F^(3,Q 2 ) , (3.22) 

where the reggeon flux f R (xp) is assumed to follow a Regge behaviour with a linear trajectory 
«k(*) = Or(0) + a' R t: 



f K M = i dt , (3-23) 

Jt cut Xp ^ ' 

where \t m i n \ is the minimum kinematically allowed value of \t\ and £ cut = —1 GeV 2 is the limit of 
the measurement. The values of B R and a' R are fixed with data from hadron-hadron collisions 
PJ. The reggeon structure function is assumed to be the pion structure function |T7|]. The free 
parameters for this component are the reggeon normalisation iV* and exponent qr. The final 
parametrization used for the fit is the sum of the three contributions detailed above (3.17, 3.18, 
I2T1 , ET22D , namely: 

F f® = _L (N in (F° (tn) + F° {in) ) + N$f£ W + m l F? (el) + N R F®) . (3.24) 
X ( ■( 1 v " '" ' 

Only one parameter (N m ) is used for the normalisation of the inelastic component since 
F ^/F^ n) = 2/9. 
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4 Phenomenology 



4.1 Fits to the HI and ZEUS data 

A fit to the recently published HI j lj and ZEUS || diffractive structure function data is 
performed separately. The result of both fits is shown in Figure 2 for HI and Figure 3 for 
ZEUS, and their parameters are given in Table 1. Let us comment the interesting features 
of the fits. The fit to the HI data leads to a very good x 2 (1-17 per degree of freedom with 
statistical errors only). The value of ckr is consistent with the usual values found for secondary 
reggeon contributions if interference effects are taken into account The value of «p is found 
to be consistent with the expected intercept for a hard BFKL pomeron fL2"| . This intercept is 
higher than the value obtained from the fit to the structure function F2 H. Qo is a typical non 
perturbative scale for the proton and very close to the value obtained in the proton structure 
function fit. It should be noted that the scale Qo appears in a quite non trivial way as the 
virtuality in the inelastic component (Q/4Qo), and in the elastic one (Q/2y r /3Qo). The validity 
of these results can be checked by starting with two different values of Qo for each component 
and the result of the fit leads to same values within errors, with the same x 2 ■ Furthermore, 
imposing a constant scale for the elastic component, i.e. Q/2y/]3Qo — > Q/Qo, leads to a very 
bad quality fit. The other features of this description will be discussed further on. 

The fit to the ZEUS data leads to a worse x 2 (x 2 /dof = 1.95 with statistical errors only). 
The result is shown in Figure 3 in continuous line, and by comparison, the dashed curve 
corresponds to the parameters obtained in the HI fit. It can be noticed that a difference is 
clearly seen both in the low Q 2 , low f3 and the high Q 2 , high {3 bins, in opposite way, the HI 
fit going from below to above the Zeus fit. 

In order to investigate the origin of these differences, a direct comparison between ZEUS 
and HI data has been performed |IB |. The HI data have been interpolated to the ZEUS closest 



bins in f3 and Q 2 using the dipole model fit. This interpolation is weakly sensitive to the model 
used as the interpolation in the kinematical variables is very small. It was checked that the 



use of the model by Bartels et al. [19[ gives a similar result. The result of this comparison is 
displayed in Figure 4. The striking feature is that the main difference between both fits noticed 
in Figure 3 comes from the region where the data are most different. The differences occur 
both at small and high f3 (i.e. small and high masses) and the solution of this puzzle cannot be 
due only to the difference in selecting diffractive events - either the rapidity gap selection [|TJ 
or the Mx substraction method @ - in both experiments. In fact, only high- mass events could 
be affected by these different selections. A global shift of ZEUS data of about 15% would be 
possible as they do not correct for proton dissociation but the difference between HI and ZEUS 
data is not only due to normalisation effects. Indeed, it is clear in Figure 4 that, depending on 
the values of Q 2 and /3, HI points lie either higher or lower than ZEUS data. The analysis of 
new accumulated data in this kinematical region would help solving the experimental puzzle 
and would be of great interest. 
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4.2 Effective pomeron intercept 



Let us now study the dependence on xp which is directly related to the rapidity gap dependence. 
One can define an effective pomeron intercept in the following way: 



.// = I (dlnFf_ l 



2 



where the t dependence is integrated out (the data are mainly at t ~ 0). This effective exponent 
can also be determinated for the inelastic and elastic components separately as a function of 
xp (Figure 5) or Q 2 (Figure 6). 

In Figure 5, The effective intercepts of both components and their sums have been compared 
to the soft pomeron intercept |2D| (dashed line) and to the bare hard BFKL pomeron obtained 
in the fit (ap = 1.395, cf. Table 1). The range of obtained values sits essentially between these 
two limits except in the large xp region (xp > 10~ 2 ). It is clearly not consistent with the soft 
pomeron value (1.08). It is as well much lower than the bare pomeron intercept || [TI| . This 
can be explained by the large logarithmic corrections induced by the a 3 (xp) term, proportional 
to log 3 (l/xp), present in both diffractive components (see formulae 3.17, 3.18, |3.21| ). The effect 



of this logarithmic term induces also an xp dependence of the intercept. Moreover, in Figure 
5, it can be seen that the xp dependence of the intercept is different between the elastic and 
the inelastic components. This induces a breaking of factorisation directly for the diffractive 
components of this model, which comes in addition to the known factorisation breaking due 
to secondary trajectories. In addition, the elastic component itself does not factorise, as it can 
be seen in Figure 5 for different Q 2 values. This comes from the mixing of the xp and Q 2 
dependence in formula 3.17, 3.18. 

As seen in Figure 6, the Q 2 dependence of the pomeron intercept (taken for instance at 
xp=0.001 where the reggeon contribution is negligible) is weak. It is in good agreement with 
the HI determination (ap = 1.204 ± .02) indicated in Figure 6. It is thus intermediate between 
the soft and hard pomeron lines but inconsistent with both of them. In our model, this softening 
of the bare pomeron comes from the large logarithmic corrections which cannot be neglected 
in perturbative calculations as was mentioned in the last paragraph. 



4.3 Scaling violations 

One striking feature of the diffractive proton structure functions was the Q 2 dependance at 
fixed xp as a function of (3 as was pointed out experimentally by the HI collaboration IJ (see 
Figure 7) and confirmed at lower Q 2 |^T| . The structure functions are increasing with Q 2 even 
at very high (3 (see Figure 7) at variance with the behaviour of the total proton structure 
function as a function of x . In the QCD dipole model, this experimental feature is described 
by a non trivial interplay between the two diffractive components. In Figure 7, the dipole fit 
is compared with the HI result showing the contribution of each component: at small /3, the 
inelastic component dominated and vary quasi linearly in log Q 2 , and at high (3, this component 
is depressed similarly to the total structure function, but is progressively substituted by the 
elastic component. Note that the enhancement present in the data at high (3 and low Q 2 which 
is probably due to vector meson production is partially reproduced by the model. It is expected 
that including the specific vector meson contributions in addition to the elastic qq component 
will even improve the result of the model in this kinematical range. It is striking that it is 



9 



possible to describe the observed scaling violations in a very different framework as the one 
given by the DGLAP evolution as was performed by the HI collaboration |T[. 



4.4 j3— dependence 

In Figure 8 is displayed the (3 dependence of the proton diffractive structure function for a 
fixed value of xp = 0.003 for different values of Q 2 . The f3 dependence is quite weak in all 
Q 2 bins and correctly reproduced by the interplay of the two components of the QCD dipole 
model. While at low Q 2 , this effect is essentially due to the inelastic component, at high Q 2 
the interplay between both components is required to describe the observed (3— dependence. 

The interplay of the different components can be analysed in more details in Figure 9 where 
the elastic, inelastic, reggeon components and their sum are displayed for three different values 
of xp and four different values of Q 2 as a function of (3. The reggeon component is only 
important at high value of xp as expected and dominates at low (3. It disappears at smaller xp. 
At low Q 2 , the inelastic component dominates in almost the full (3 range, while at higher Q 2 , 
it is only important at low (3 and replaced at low masses by the inelastic component. 



4.5 Longitudinal contribution to the diffractive structure function 

In Figure 10 are displayed as a function of (3 in different Q 2 , xp bins the longitudinal and 
transverse components of the proton diffractive structure functions. Note that the xp values 
are chosen a bit smaller than in the previous figure to depress the reggeon contribution. Note 
also that we did not separate the longitudinal and the transverse contributions for the inelastic 
component as they are directly proportional by a factor H L {l/2) / H T {l/2) = 2/9 (see formula 
|3.21|) . In addition to the already mentioned dominance of the inelastic component at small 



Q 2 , the longitudinal elastic component is found to be high at high (3 and crosses over the 
transverse component near (3 ~ 0.8. It is thus expected to obtain high values of the ratio R 
of the longitudinal to the transverse components at high values of (3. The result of our fit 
prediction is displayed in Figure 11. We note that the R ratio remains small (~ 0.2) in almost 
the full kinematical plane except notably at high (3 where it may reach high values such as 2. 
Note that this value is in the range of the measured R ratio with vector meson production || . 
A measurement of R in diffraction would thus be of great interest and would be a good test 
of the model. The question arises whether such high values of R would not modify the _F 2 D 
measurement itself as it was assumed to be HJ . We have checked that the parameters and the 
quality of the fit are not much influenced if only data points with y < 0.3, where the influence 
of R is small, are taken. 

It is instructive to notice that another model of diffraction based on selecting qq and qqg 
components of the photon |L9] also leads to a large contribution of the longitudinal qq contribu- 



tion at high (3. However, a distinct feature of the dipole model is the above mentioned difference 
between the high and low Q 2 behaviour where at low Q 2 , the inelastic component dominates 
and induces a small value of R. We thus expect thet a measurement of R in diffraction at low 
Q 2 would be a way to distinguish both models. 
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5 Conclusion and outlook 



Let us summarize the results we have obtained. First of all, the colour dipole model has provided 
us with a interesting theoretical and phenomenological framework to study hard diffraction. On 
thetheoretical ground, it enables to interpret hard diffraction in terms of BFKL dynamics of 
(resummed) perturbative QCD, which allows one to make theoretical estimates. On the other 
hand, the success of the phenomenological fit using analytical amplitudes allow a discussion of 
the various interesting and intringuing aspects of the data. 

The model predicts two contributions: an elastic one for which the final state consists in 
the hadronization of a quark-antiquark pair, and an inelastic one, including soft gluon radia- 
tion reproduced in terms of a cascading process involving colour dipoles. Both contribute in 
quite different kinematical domains since they are characterized by different mass distributions 
(different f3) and for different photon polarization. In this paper, we have obtained a suitable 
parametrization and a satisfactory fit to the HI data along these lines. The fit to Zeus data is 
less satisfactory and with different parameter values. By a model-independent analysis of both 
data, we could identify the reason of the discrepancy in a genuine difference between data sets 
when interpolating in the same bins. This point deserves a detailed experimental study in the 
near future, since it may be of importance to doiscriminate between models, in particular the 
confirmation of the trend of present HI data would be in favor of the dipole model. 

On the theoretical ground, we have computed amplitudes obeying some of the requirements 
of BFKL dynamics for diffractive amplitudes (on a dipole). However, exact calculations of these 
rather complicated amplitudes are still under way [^] and we hope in the near future obtain 
definite predictions with all requirements of BFKL dynamics, the evaluation of next-leading 
BFKL effects would also be welcome since they already play a role for the total structure 
functions. 



Finally, as emphasized in ref. |TT[ the colour dipole model formalism calls for a unified 
description of the diffractive and total deep-inelastic scattering events, e.g. including events 
with no rapidity gap. We showed that within the precision of the current data, there are quite 
a few indications (similar scale Qq, softening of the gard Pomeron by logarithmic facors in 
diffraction,etc.) of such a common theoretical ground. However, further tests of the model 
are deserving. The first one would be a confrontation of the predicted R ratio with the data 
if available: indeed, the various models should predict quite different contributions from the 
two polarization states of the photon. Other useful tests concern the final states. For instance 
one can compute the predictions for diffractive vector meson production and confront them to 
the recent data. Such tests might help distinguish between the few different models for hard 
diffraction which are able to describe the data, definitely, the QCD dipole model of BFKL 
dynamics is one of them. 
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TABLE CAPTION 



Table I 

Parameters obtained for the fit. 

The fit has been performed with statistical errors only. The first error quoted is the statis- 
tical one and the second the systematic one. 





HI 


ZEUS 


ap 


1.395 ± 0.005 ± 0.003 


1.327 ± 0.001 ± 0.016 




0.682 ± 0.046 ± 0.049 




Qo 


0.428 ± 0.011 ± 0.001 


0.241 ± 0.002 ± 0.014 




0.00244 ± 0.00028 ± 0.00033 


0.00374 ± 0.00015 ± 0.00076 


Nj! 


40.0 ± 1.6 ± 3.6 


126.9 ± 0.5 ± 37.6 


Nf 


12.8 ± 1.2 ± 1.1 


24.43 ± 0.51 ± 8.12 


N R 


7.44 ± 1.84 ± 3.50 


± 58 ± 


x 2 


255.4 for 226 pts 


89.6 for 53 pts 




1.17/d.o.f. 


1.95/d.o.f. 



Table I 
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FIGURE CAPTIONS 



Figure 1 Schematic description of the QCD dipole model of diffraction. 1-a: elastic component, 
1-b: inelastic component 



Figure 2 Result of the F® fit to the HI data j^]. The data are displayed by triangles (with 
statistical and systematic errors added in quadrature) as a function of xp in (3 and Q 2 bins. 
The fit has been performed with statistical errors only and is displayed in full line (see text). 



Figure 3 Result of the fit to the ZEUS data ||/. The data are displayed by triangles with 
statistical and systematic errors added in quadrature as a function of xp in (3 and Q 2 bins. The 
fit has been performed with statistical errors only and is displayed in full line (see text). For 
comparison, he result of the fit of the HI F® data is displayed in dashed line. 



Figure 4 Direct comparison between HI (squares) and ZEUS (triangles) data. The HI data 
have been interpolated to the ZEUS bins by using the dipole parametrisation. Both data 
are found compatible within error bars. Note however that sensitive differences are seen in 
particular in three different (3 and Q 2 bins (Q 2 =8 GeV 2 for (3=0.2, Q 2 =Q0 GeV 2 , for (3=0.7 
and 0.9) (see text). 



Figure 5 Effective pomeron intercept as a function of xp. The dependence of the effective 
pomeron intercept ctp on xp is shown for the inelastic component (independent of Q 2 ; continuous 
line) and for the inelastic one for 3 different values of Q 2 . The value of (3 is fixed at 0.3. The 
straight lines at ap=1.395 and ap=1.08 correspond respectively to the bare BFKL intercept as 
found by the fit (see table 1), and the soft Donnachie Landshoff prediction [EH . 



Figure 6 Effective pomeron intercept as a function of Q 2 . The effective pomeron intercept of 
the model is displayed: elastic component (dashed curve), inelastic component (dotted curve), 
total (full curve). The value of (3 and xp are fixed for reference at ((3 = 0.3 and xp = 10~ 3 ). 
The bare and soft pomeron intercepts are given by dotted-dashed straight lines. The HI 
determination [[U lies in the region between the two dashed straight lines. 



Figure 7 Scaling violations. The dependence of xpF 2 D on Q 2 for different values of (3 at fixed 
xp (3.10 -3 ) is shown together with the dipole model fit. Dotted line: elastic component, dashed 
line: inelastic component, full line: total. 
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Figure 8 Beta dependence. The dependence of xpF 2 D on (3 for different values of Q 2 at fixed xp 
(3.1CT 3 ) is shown together with the dipole model fit. Dotted line: elastic component, dashed 
line: inelastic component, full line: total. 



Figure 9 Different components of the HI F 2 D fit. Dotted line: elastic component, dashed line: 
inelastic component, dashed-dotted line: secondary reggeon component, full line: sum of all 
components. 



Figure 10 Different components of the HI F® fit. Longitudinal versus transverse. Dotted line: 
tranverse elastic component, dotted-dashed line: longitudinal elastic component, dashed line: 
total (Ft + F L with F L /F T = 2/9, see text) inelastic component, full line: result of the fit. 



Figure 11 Prediction for R, the ratio of the longitudinal to the transverse cross section. Note 
the different scale for the high /3 bins. 
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